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ABSTRACT Close to melting transitions it is possible to propagate solitary electromechanical pulses which reflect 
many of the experimental features of the nerve pulse including mechanical dislocations and reversible heat production. 
Here we show that one also obtains the possibility of periodic pulse generation when the boundary condition for the 
nerve is the conservation of the overall length of the nerve. This condition generates an undershoot beneath the baseline 
('hyperpolarization') and a 'refractory period', i.e., a minimum distance between pulses. In this paper, we outline the 
theory for periodic solutions to the wave equation and compare these results to action potentials from the femoral nerve 
of the locust (locusta migratoria). In particular, we describe the frequently occurring minimum-distance doublet pulses 
seen in these neurons and compare them to the periodic pulse solutions. 
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Introduction 

In several recent publications we have proposed that the ac- 
tion potential in nerves is an electromechanical solitary wave 
or soliton |Q] [2] [3] |4) . Although the propagation of mechanical 
pulses in nerves has been discussed before 16] 13 [H |9), we 
have provided an quantitative formalism based on thermody- 
namics for how such pulses are generated and how they propa- 
gate. An important justification for the assumption of an elec- 
tromechanical process is the experimental observation of re- 
versible heat changes in phase with the action potential and a 
zero net heat release during the action potential (e.g., ifTOlfTTl 
|T2j |T3] [T4) ) . This finding is incompatible with an electrical pic- 
ture of the nerve pulse (Hodgkin-Huxley model [ 15 1) in which 
currents (of ions) passing through resistors (i.e., ion channel 
proteins) play a dominant role. The latter model would result 
in a net heat release that is not found in experiments. The ini- 
tial heat release during the action potential was also found to be 
much larger than the free energy required to charge the mem- 
brane capacitor, which rules out the possibility that a purely 
electrical picture is able to explain the nerve signal. The pic- 
ture of the nerve pulse as an electromechanical phenomenon 
is further supported by the finding of mechanical changes of 
the nerve under the influence of the action potential, including 
changes in thickness and length (e.g., llr^ [T6l[T4ll ). 

The dynamic equation underlying the soliton description is 
based on the simple wave equation for sound with the lateral 



density of the nerve membranes as a variable. When describing 
sound propagation in air, it is usually assumed that the sound 
velocity, Co, is constant. This leads to the simple differential 
equation d 2 Ap/dt 2 = Cq<9 2 Ap / 'dx 2 . Close to the melting 
transition in membranes, however, the speed of sound is a 
sensitive function of density ifTTl [T8l [T9l |20 1 and frequency 
ETlfTTIl . The simultaneous presence of non-linearity and dis- 
persion gives rise to the possibility of localized excitations or 
solitons. The resulting differential equation is somewhat more 
sophisticated but now displays analytical localized solutions 
E2l |2) Q Due to the fact that biomembranes carry charges, 
these solitons are of an electromechanical or piezoelectric 
nature. Thus, they are not only localized voltage pulses but 
also display reversible mechanical features and a reversible 
change in heat that are in agreement with experiment. Further, 
since our theory is of a thermodynamic nature, it implicitly 
has a role for all thermodynamic variables, e.g., temperature, 
pressure, electric potential and the chemical potentials of drugs 
and ions. This implies that the generation of the nerve pulse 
can be affected by pH 1231 . calcium concentration, and general 
or local anesthetics f2l l23"l . 

While a reversible electromechanical picture of the nerve pulse 
is compelling and able to address many more physical features 
of the nervous impulse than traditional models, there remain 
issues that have not yet been addressed. In particular, this in- 
cludes refractory periods, hyperpolarization, and pulse trains. 
The refractory period is the minimum possible time between 
two nerve pulses. In the Hodgkin-Huxley picture it is a conse- 
quence of the complicated voltage and time dependence of the 

Here, we use the term 'soliton' as synonymous with 'solitary wave'. Gen- 
uine solitons can pass through each other without dissipation, which is not the 
case for our pulses [22]. Strictly speaking, the term solitary wave is thus more 
appropriate. 
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opening characteristics of ion channels. In simplified terms, the 
refractory period is due to the relaxation time of the ion chan- 
nel proteins after firing and is several milliseconds long, i.e., 
several pulse widths. Hyperpolarization is the undershoot of 
the voltage below the resting potential during one phase of the 
action potential. 

In the soliton model the action potential consists of a local re- 
gion of the nerve with high chain order and high local den- 
sity. An immediate consequence is that the excited nerve con- 
tracts. This has been observed in isolated nerves |6][14). How- 
ever, it is difficult to imagine that nerves also contract in the 
native situation when they are embedded in tissue where dis- 
tances are kept constant. Here, we demonstrate that one can 
obtain solutions for the electromechanical differential equation 
that maintain the length of the nerve. They consist of periodic 
pulses with an intrinsic minimal distance between pulses that is 
of the order of about 5-10 pulse widths. Using locust femoral 
nerves, we demonstrate that one finds minimal distances of ex- 
actly that order. In particular, we describe a phenomenon that 
we called pulse "doublets" that displays a constant distance be- 
tween pulses even if measured under different conditions. We 
have found such doublets in locust femoral nerves and also 
recently doublets and triplets in crayfish motor neurons. We 
show that the refractory period in the soliton model is a con- 
sequence of the conservation of mass and overall length of the 
nerve. This boundary condition also leads to periodic solutions 
or pulse trains and to hyperpolarization. 

Materials and Methods 

Locust nerve recordings: Adult male and female Locusta 
migratoria were reared in a laboratory colony at about 30° C 
and with a 12:12 hours dark-light cycle. Experiments were 
performed using amputated metathoracic legs which were 
mounted on a platform and fixed on its external lateral face 
with PlasticineTM (Modeling Clay KR 25, non-toxic, Becher 
Gottingen). The axis of rotation of the femoro-tibial joint was 
aligned with the origin, of a semicircular angular plane. Inside 
a Faraday cage and with the aid of an optical microscope, a 
small window of the postero-lateral cuticle was removed in 
order to expose the peripheral lateral internal nerve (properly 
called posterior tegumental femoral nerve n5B2c l24ll25l . and 
the tibia was cut to a stump to allow oxygenation inside the leg 
via the trachea. See figure [3] (bottom). In order to record the 
electrical activity from the single antero-lateral (RDAL) and 
the two postero-lateral multipolar sensory afferents (RDPL and 
RVPL) (nomenclature according to [26]), monopolar platinum 
wire electrodes, 50 /im in diameter, were hooked on the lateral 
internal nerve n5b2c. 

With the aid of a pipette, a small drop of locust saline 
(NAC1 140 mM, KC1 10 mM, Na H 2 P0 4 5mM, CaCl 2 
2 mM,Saccharose 90 mM; pH 6.8) was deposited on the 
wounds to avoid the coagulation of the hemolymphe and also 



used as a reference. A small amount of Vaseline (DAB 10, 
Roth GmbH, Karlsruhe, Germany) was added to isolate the 
nerve-electrode system as well as to avoid desiccation of the 
preparation. 

The femoro-tibial angle was extended and fixed with an ento- 
mological pin taking care to avoid lateral strain on the joint. 
Although primitive, the use of a pin to fix the angle is the most 
precise method, because the tibia tends always to the flexion 
and this counterforce allows a perfect static equilibrium at the 
desired angle of extension. All experiments were performed at 
room temperature of 22-24° C. 

Detailed Procedure: 

Results presented in this report belong to 10 samples taken 
from different animals. In all cases, three electrodes were 
hooked on the lateral internal nerve n5b2c, and a small amount 
of Vaseline was added to isolate each electrode-nerve subsys- 
tem. The distance from the axis of rotation of the joint, O, to 
the electrodes was 10, 12 and 15 mm, respectively. We did 
not put the first electrode closer to the origin of the system for 
fear of damaging the receptors, and the last electrode could 
not be hooked beyond 17 mm because the nerve under study 
joins the nerve n5b2 at that point. In all experiments, the tibia 
was extended from to 80° to verify that no activity could be 
recorded. Afterwards, at each angle from 80 to 150° in steps of 
10°, the tibia was fixed completely with an entomological pin 
and sessions of three minutes were recorded for each angle. 
We found that 150° is the limit of extension since beyond this 
angle the end of the tibia causes pressure against the distal 
dorsal edge of the femur and therefore the conditions are no 
longer physiologic [26|. All experiments were performed at 
room temperature of 22-24° C. 

The action potentials were picked up externally by platinum 
wire hook electrodes. A preamplifier (four channels, lOOOx 
amplification, no pass filter, 9V battery source, workshop of 
the Institute for Zoology, University of Gottingen) was used 
for primary signal amplification within the Faraday cage, and 
secondary amplifier (four channels, lOx amplification, 100 Hz 
high pass filter, 30 kH low pass filter, 15 V transformer source) 
used to pre-filter the signals. The amplified action potentials 
were digitized via PC based data acquisition hardware (Packard 
Bell desktop computer, 8 channel RUN Technologies acquisi- 
tion port, PCI-DAS 1200Jr board, Microsoft Windows XP Pro- 
fessional SP2, 2.66GHz Pentium(R) 4 CPU, 512 MB RAM) 
and recorded and analyzed with the DATAPAC 2K2 software 
(RUN Technologies, Mission Viejo, CA, USA). 

Analysis of the recordings: 

The smallest digital recording sample period achieved on the 
system was 49.6 /is; periods of 89.6 and 148.8 /is were cho- 
sen for longer recordings. No software-based secondary filter- 
ing was performed on the recorded channels. Pulse height was 
defined by the voltage distance between maximum depolariza- 
tion and maximum hyperpolarization in order to account for a 
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potential signal baseline slope during the recording. Pulse du- 
ration was defined as the time difference between maximum 
depolarization and maximum hyperpolarization during a pulse. 
In all pulse trains, the onset of the individual signals was de- 
fined by the maximum depolarization within the pulse, the off- 
set by the maximum hyperpolarization. Signal velocity along 
the nerve was estimated by subtracting the onset times of time- 
linked signals between the channels and multiplying the inverse 
by the measured distance between the corresponding hook elec- 
trodes. The instantaneous pulse period was defined as the pe- 
riod between the current onset and the onset of the next signal 
within a pulse train. 

In each region of interest, pulses (APs) were first identified and 
time localized via a simple threshold pass analysis. The thresh- 
old was established visually, positioned at roughly 80% of the 
maximum depolarization amplitude of the smallest unit of in- 
terest, then adjusted manually in such a way as to maximize 
the sensitivity and specificity regarding AP localization, ideally 
with no noise spikes falsely identified as APs. The onset of the 
marked pulses was then moved to the first peak depolarization, 
the offset to the first peak hyperpolarization. 



Theory 

In the following we want to derive the equations of motion in 
a cylindrical membrane along the coordinate x. The lateral 
density, p A , of a membrane is larger in the gel phase of the 
membrane than in the fluid phase because the area per lipid 
molecule changes by about 24%. For this reason, one can in- 
duce a phase transition not only by adjusting the temperature 
but also by changes in lateral pressure or other thermodynamic 
variables. 

Let us consider the nerve axon as a one-dimensional cylinder 
with lateral density excitations moving along the coordinate x. 
As outlined in (TJ, sound propagation in the absence of disper- 
sion is governed by the equation 
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where Ap = p — p A is the change in the area density of 

the membrane as a function of x and t and c = \J\/(n A p A ) is 

the density-dependent velocity of sound. Here, p A is the den- 
sity of the membrane at physiological conditions and is slightly 
above the melting transition. The above equation is related to 
the Euler equation in hydrodynamics and will not be derived 
here. To the extent that the compressibility is independent of 
the density and the amplitude of the propagating density wave, 
Ap A << p A , the sound velocity is approximately constant 
(c = cq). Eq.Q reduces to 
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This is the well-known equation that governs sound propaga- 
tion in air. However, in several publications we showed that the 
compressibility, k a , depends sensitively on temperature (and 
therefore also on density) close to the melting transition of lipid 
membranes ifTTl [T8l [T9l l20l . Therefore, the above simplifica- 
tion cannot be made. Both the liquid and gel phases are rel- 
atively incompressible. At densities near the phase transition 
where the two phases co-exist, a small increase in pressure can 
cause a significant increase in density by converting liquid to 
gel. Near this phase transition, the compression modulus is 
dramatically smaller. We thus approximate the sound velocity, 
c, as 

1 
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with p < and q > 0. The sound velocity is also frequency 
dependent [21 ]. This means that the system is dispersive, which 
is a necessary requirement for soliton formation. At very low 
frequencies, lu 0, the adiabatic compressibility approaches 
the isothermal compressibility. 

In the following, we use the values for the isothermal com- 
pressibility as a low frequency approximation of the the adia- 
batic compressibility. For unilamellar DPPC vesicles we found 
c = 176.6 m/s, p = -16.6 c 2 / p A and q = 79.5 c 2 ,/ (p A ) 2 
with p A = 4.035 x 10 -3 g/m 2 , assuming a bulk temperature of 
T=45°C (TJ. Similar values are obtained for biological mem- 
branes at physiological temperature. 

The isentropic compressibility decreases at higher frequencies 
and leads to a propagation velocity which increases with in- 
creasing frequency. We approximate such dispersive effects by 
introducing a term, —hd 4 Ap A /dx 4 with h > 0, in eq. 111. As 
demonstrated in [1 1, this term implies a linear dependence of 
the pulse propagation velocity on frequency. This form of the 
dispersion term represents the leading term in a typical long- 
wavelength dispersion relation. In the absence of quantitative 
knowledge of the dispersion present in nerves, this term has the 
advantage of being relatively passive in that changes in h alter 
the spatial size of solitary waves but not their functional form. 
The inclusion of additional terms can lead to soliton instability 
and should not be considered without empirical justification for 
the form chosen. 

The equation governing sound propagation now becomes 
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This equation is closely related to the Boussinesq equation 
[27 1. Since we are interested in finding solutions that propa- 
gate without changing their shape, we transform to a moving 
coordinate system with z = x ~ vt. The time- and space- 
dependent differential equation (eq. [4]) can be transformed into 
a time-independent equation by using dAp/dx = (dAp/dz) ■ 
(dAz/dx) = dAp/dz and dAp/dt = (dAp/dz) ■ (dz/dt) = 
—vdAp/dz. 
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Figure 1: The potential V(Ap A ) from eq. (|s])/or three different integration constants, C, and a velocity v — 0.8 • Co- The 
inserts show the potentials in a larger density range showing that they all possess two maxima and one minimum. The 
electromechanical pulses correspond to a movement in the potential between the limits given by the two dashed lines. Left: 
For C = g/cm 2 one finds a singular solution given in Fig.\2\ right top. This solution displays a baseline of Ap A = 0. 
Center: For C — 2 g/s 2 one obtains a periodic solution Fig.\2\ right, 3rd trace from top. This solution displays a baseline 
ofAp A < 0. Right: For C = h.!2g/s 2 one obtains a potential with two maxima of equal height, corresponding to the 
periodic solution given in Fig^^ right bottom. This solution also displays a baseline of Ap A < 0, and the peak hits a 
maximum and broadens. 



Eq.Q can now be rewritten as 
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We now search for periodic solutions. We can now perform two 
integrals with respect to z from one minimum of the periodic 
solution to the next 
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This can be easily checked, since the second derivative of 
eq. (|6]l yields eq. ([5j. C is an integration constant. The second 
integration constant disappears for reasons of symmetry. 

Now multiply both sides of eq. ^ by d(Ap A ) / dz and integrate 
once more to yield 
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Here, Vq is a further integration constant. In this equation there 
is one term that is proportional to {dAp A /dz) 2 and a second 
term that explicitly only depends on Ap A . Since a mechanical 
analogy is useful, we call the first term a "kinetic term" and the 
second one a "potential term": 



A/ 



(8) 



'kinetic term' 



This equation is formally equivalent to the equation of mo- 
tion \mv 2 + V(x) = const, familiar from classical mechan- 
ics where z plays the role of time and the density change Ap A 
serves as a spatial coordinate. The shape of the soliton is a re- 
sult of the motion in this potential, and all solutions of eq. d9]l 
can be classified according to this potential energy surface [28 ]. 
It should be stressed that this analogy to classical mechanics 
must not be taken too literally. It is just an analogy because 
one of the terms depends on the square of the derivative of the 
density with respect to the coordinate, while the second term 
is an analytic function of the density. Vo is not relevant be- 
cause it moves the potential up or down without any influence 
on the solutions of the differential equation. The value for the 
integration constant C will depend on experimental boundary 
conditions. 

In the following we explore the solutions of eq ([9]). 



Results 

Solitary and periodic pulse solutions: 

Boundary conditions are needed to determine the desired solu- 
tions. In the following we assume that the total length of the 
nerve is fixed. This implies that the mean area density of the 
nerve membranes stays constant (i.e., (Ap' 4 ) = 0), and that 
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Figure 2: Solutions to eq. ^ for different values of the velocity vand the integration constant C. Left: For C = Og/s 2 one 
obtains solitary pulses. Their peak amplitude depends on the velocity, v. Right: Solutions for v = 0.8 • Co and various 
values of the integration constant C. For all C > Og/s 2 the solutions are periodic and one finds an undershoot. The pulse 
distance is a function of C. 



positive peaks must coexist with regions of negative density 
change. As shown below, this condition leads to periodic pulse 
solutions with an undershoot of the pulse (hyperpolarization) 
and refractory periods. We use the elastic parameters for the 
lipid DPPC given above. However, as shown in HI values 
appropriate for biological membranes would work quite as 
well. The value for C will be varied. For instance, for a fixed 
propagation velocity v = 0.8co one obtains stable solutions in 
the range g/s 2 < C < 5.72 g/s 2 . 

The potential V(Ap A ) is a polynomial of order 4 with negative 
sign in the highest order term. Therefore, this potential can dis- 
play two maxima and one minimum. Stable solutions are con- 
fined to be between the two maxima; all other solutions would 
be unstable. For C = g/s 2 the potential is shown in Fig.[T[ 
left. It displays one local maximum at Ap A = and a second 
maximum at positive Ap A with higher value. The shape of the 
soliton is defined by the movement within the potential indi- 
cated by the arrow. It is defined by moving within the potential 
from the local maximum at Ap = to that positive value 
for Ap A for which the potential assumes the identical value. 
This defines the amplitude of the soliton. One obtains the iso- 
lated solitary solutions shown in Fig. [2] left. For each ampli- 
tude one obtains exactly one soliton with well-defined velocity 
somewhat smaller than the small-amplitude speed of sound, Co. 
Larger amplitude results in slower propagation velocity. There 



is a limiting velocity, vn m i t = c 2 — p 2 /6q ~ 0.65 • cq at which 
the solitons assume maximum amplitude and minimum veloc- 
ity HI. The solutions consist of a local change in density above 
the baseline. The neuron is considered being infinitely long. 
Therefore the mean density is still the baseline value even if it 
is locally different from zero. The solitary character is a con- 
sequence of the potential having zero slope at Ap A = 0. If the 
movement of the membrane starts from a value of Ap A > 
one obtains periodic solutions. However, those solutions dis- 
play a mean density change larger than zero and are therefore 
no solutions are possible under the assumptions of constant 
nerve length. 

For < C < 5.72 g/s 2 , the first maximum of Vis located at 
a negative value of Ap A . This is shown in Fig.[T] center. Un- 
der these conditions no solitary solution can exist because that 
would result in (Ap A ) < 0. However, there are periodic so- 
lutions that display a mean density change of zero. Such solu- 
tions are shown in Fig.|2](right) a fixed velocity of v — 0.8 ■ cq. 
For each value of C and given velocity v pulse separation and 
amplitude are uniquely defined. 

The upper limit for the integration constant is C « 5.72 g/s 2 . 
This case is shown in Fig. [T] (right). One recognizes that for 
this value of C the two maxima of the potential assume the 
same value. This influences the shape of the pulse such that it 
becomes flat on the top (see Fig. [2] (right), bottom trace). No 



5 



stable solutions exist for C > 5.72 g/s 2 and for C < Og/s 2 . 
In summary, one finds the following: 

1 . The integration constant C depends on the experimental 
conditions. There exist solitary and periodic solutions for 
the differential equation [9] and the boundary condition of 
constant mean density. 

2. For each C and each velocity v there is exactly one solu- 
tion with a mean density change of zero. For C — g/s 2 
this is a singular pulse (soliton) and for C > g/s 2 it is 
a periodic pulse. 

3. There is a minimum distance between pulses dictated by 
the requirement that the mean density is constant, corre- 
sponding to a refractory period (Fig. [2). 

4. For each pair of pulse amplitudes and pulse distances 
there is exactly one velocity. 

This implies that the velocity of the pulse train depends on am- 
plitude and frequency of the stimulation process. It further im- 
plies that any change of the amplitude resulting from a change 
in the conditions (e.g., the presence of anesthetics, changes in 
temperature or pressure, etc.) will have a predictable conse- 
quence for the propagation velocity and the spatial distance 
between pulses. This means that the individual pulses of the 
periodic solutions influence each other. One further result from 
Fig.|2]is that increasing pulse amplitude first leads to a narrow- 
ing of the pulses. When the pulses approach their maximum 
amplitude, they become broader and display a plateau at max- 
imum density. Pulse shapes of this kind are known in nerves, 
for example the prolonged action potentials of squid axons that 
occur under perfusion with aqueous media containing unusual 
salt composition |29l . 

The above considerations apply to an axon that is infinitely 
long. In such a system, there are no constant velocity, shape- 
preserving solutions involving a finite number of pulse peaks. 
Solutions in which a finite segment of a periodic solution is 
joined smoothly to the solution Ap A = is expected to be 
long-lived since only the limited material at the interface can 
lead to its disintegration. Of course, real nerves have finite 
length and only a limited number of pulses can fit into a neuron 
at a time. It is therefore to be expected that the situation in the 
real nerve is close but not identical to the idealized situation 
described above. In particular, we suggest that the pulse dou- 
blets observed in the locust femoral nerve (described below) 
represent a limiting case of the periodic solutions considered 
here. The existence of periodic solutions, such as those shown 
in Fig.|2j is a very robust result. Arguments similar to those 
made here can be adopted to find periodic solutions for essen- 
tially any system exhibiting isolated solitons. This includes, 
e.g., the FitzHugh-Nagumo model 11301 [3D which is closely 
related to the Hodgkin-Huxley model. The existence of (ap- 
proximate) periodic pulses is thus not likely to be useful in 
discriminating between models of the action potential. In this 




Figure 3: Representation of the locust metathoracic leg 
showing the main nerves and specifically the lateral 
nerve n5b2 we studied (image adapted from 138]). The 
femoral nerve extends along the whole leg. The angle of 
leg extension, a, influences firing frequencies ( cfi, Fig. |6|. 



connection, we note that the model employed here was devel- 
oped to describe myelinated nerves. In HI we suggested that 
myelination restricts the nerve to longitudinal density perturba- 
tions. This allowed us to adopt the one-dimensional represen- 
tation of the propagating pulse given in eq. [8] above. Without 
a myelin sheet we expect that the mechanical dislocations will 
contain out-of -plane components which would render the pulse 
significantly slower. So far, we have no mathematical model 
for this scenario. However, we expect such pulses to be based 
on physics similar to the longitudinal pulses described above. 
Below, we compare our theoretical results to pulse doublets in 
locust femoral nerves. 

Locust femoral nerve: Pulse trains and minimum 
distance doublets 

The femoral nerve of the locust contains three neurons (Pig.[3j 
from the extension of the femoro-tibial joint. The three neu- 
rons typically display periodic patterns with frequencies and 
amplitudes specific for each unit. The firing frequency also de- 
pends on the angle of extension, a (see Fig.[3j. Therefore, in 
a recording one can distinguish the firing patterns of the three 
neurons and analyze them separately (Fig. [4]). We have labeled 
them by signal amplitude and refer to them as neuron 1 , 2 and 
3. We have measured the pulse velocity by recording the action 
potentials at two positions of the femoral nerve separated by 
5mm (Fig. [4]). We found 

1. v = 1.02 ± 0.01 m/s for neuron 1 (large amplitude) 

2. v = 0.79±0.01m/s for neuron 2 (medium amplitude) 

3. v = 0.73 ± 0.01 m/s for neuron 3 (small amplitude) 

Most of the pulses in the three neurons have time intervals of 
several times 10 ms, depending on the angle of extension of the 
leg. Fig. [5] shows recordings that display pulse doublets with 
identical amplitudes and separations. They occur in all three 
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Figure 4: Recordings from the femoral nerve of the lo- 
cust: Three different pulse originating from three differ- 
ent neurons that are distinctly different in amplitude and 
shape. Indicated are the recordings from two electrodes 
that are separated by 5 mm. The time lags between the 
two electrodes yield the pulse propagation velocity which 
is of order 1 m/sfor all three neurons (see text). 



neurons. These pulse doublets occurred quite frequently and 
account for some 10-20% of all action potentials. The total 
length of the locust femoral nerve is about 2 cm. The pulse 
width is about 2 ms. The minimum distance of pulses we have 
recorded is about 1 1 ±2 ms (see below). In 1 1 ms the nerve 
pulse travels approximately 1 cm. This implies that one can- 
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Figure 5: Pulse doublets occurring during the record- 
ings of the femoral nerve. Top: Recording from neuron 
2. Center: Recording from neurons 1 and 2. Bottom: The 
two doublets marked in the top traces given in higher res- 
olution. 



not fit more than 2 or 3 pulses in the same neuron at the same 
time. All pulses with distances larger than 20-30 ms must be 
regarded as uncorrelated, i.e., about 80-90% of all pulses. By 
uncorrelated we mean that the pulse shapes of those pulses can- 
not influence each other because they do not coexist in a nerve. 
However, the minimum distance pulse doublets are simultane- 
ously present in the neurons and can therefore display a mutual 
influence that may influence the separation and shape of the 
pulses. For this reason we will in the following consider the 
doublets as an approximation to the periodic solutions in the 
grasshopper nerve. It should be noted that we recently recorded 
similar pulse doublets and even triplets in the crayfish motor 
neuron 3 (data will be published elsewhere). We have recorded 
a histogram of pulse intervals of neuron 2 as a function of the 
angles of extension (Fig|6]l. One finds that for each of the three 
angles, 90°, 120° and 150°, the histogram displays 2 maxima. 
We fitted them to Gaussians (Fig. [6]) and determined the mean 
peak positions and the percentage of doublet intervals (shown 
in Table[T|. About 10% of all intervals in the recordings are 
related to doublets, corresponding to about 20% of all action 
potentials being doublet peaks. While we find that the longer 
interval is strongly affected by the angle of leg extension, the 
short interval corresponding to the doublet peak intervals is un- 
affected by the angle of extension within error. For this reason 
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Figure 6: Top: Histograms of pulse separation in neu- 
ron 2 for three different angles of extension: a =150°, 
120° and 90°. Each histogram contains two peaks, one of 
which is dependent and the other one being independent 
of the angle of extension, a. Bottom: The small interval 
peak corresponds to a doublet spacing of about 11ms. 
Details are given in table^ 
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they are listed independently in table [T] The short interval is 
on average around 11 ms corresponding to the doublet spacing. 
The width/distance ratio for the doublet peaks is 0.1 — 0.2 de- 
pending on the somewhat arbitrary criterion for defining width. 
The closest distance of the theoretical periodic pulses (Fig. [2} 
corresponds to a width/distance ratio of about 0. 125, i.e., it is of 
same order. Therefore we suggest that the doublets correspond 
to the closest possible pulses in our theory. 



Table 1: Mean intervals of action potentials in neuron 2 ex- 
tracted from the peaks of the distributions in Fig. [6] 



angle 


peak 1 
position width 
[ms] [ms] 


peak 2 
position width 
[ms] [ms] 


% doublet 
intervals 
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39.7 
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82.5 
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25.8 
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all 

angles 
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Discussion 

The Hodgkin-Huxley model for the nerve pulse describes 
the nerve impulse as a voltage pulse generated by time- and 
voltage-dependent opening of ion channel proteins. As we 
have discussed before, this model is unable to explain a 
number of features of the nerve pulse including the genera- 
tion of reversible heat and the observed thickness and length 
changes H] [2] [3j O. Reversible heat is a particular challenge 
to the HH-model because it implies that the physics of the 
nervous impulse is based on thermodynamically reversible 
processes while the HH-model is exclusively based on irre- 
versible processes due its reliance on the flow of ions along 
concentration gradients. In order to address these problems 
we have proposed that the nerve pulse instead consists of an 
electromechanical density pulse (soliton) that does not rely 
on ion channel proteins. In particular, we have shown that 
close to melting transitions in biomembranes the possibility of 
soliton propagation, i.e. of localized pulses, exists. This model 
automatically addresses the observed changes in length [6| and 
IH and thickness EH [16] [33] El El ESI and the reversible 
heat production |[T0lfTT1fT3l[T4ll36l. 

The key feature of our theory is a localized pulse with con- 
stant entropy, i.e., without net heat exchange with the environ- 
ment. Since our theory is of thermodynamic nature, any change 
in a thermodynamic variable that has the potential to move 
the membrane through its transition is also able to generate a 
pulse. This includes pulse generation by changes in tempera- 
ture, voltage, lateral pressure, pH, calcium concentration, etc. 
In particular, we have shown that this theory implies a mech- 



anism to explain anesthesia [2, 23]. So far, however, our the- 
ory has not included some of the features observed for the ner- 
vous impulse including hyperpolarization, refractory periods, 
and pulse trains. These features are addressed in the present 
paper. Our original approach to pulse propagation refers to 
an infinitely long axon that is quasi-one-dimensional, mean- 
ing that the thickness of the nerve is negligible compared to 
the length of the pulse, and that the density perturbation is of 
longitudinal nature. The pulse consists of a locally compressed 
region in which the membrane is transiently pushed through 
its melting transition while releasing and reabsorbing the latent 
heat of the transition. For an infinitely long nerve the mean 
area density in the presence of a singular pulse is unaffected. 
However, for a real nerve of finite length the overall length of 
the nerve should decrease while a nerve pulse is traveling along 
the axon. Such length changes have been observes in isolated 
nerves ||6] Q3) . However, it seems unlikely that a nerve within 
a body can undergo significant length changes because they are 
attached to muscles or other nerves with more or less fixed po- 
sitions. Therefore, we must distinguish between nerves with 
constant length and nerves removed from bodies without any 
external constraint on length. 

Here, we have shown that the constraint of fixed overall length 
leads to the emergence of periodic solutions. These solutions 
share the property that their average area density is constant, 
leading to an undershoot under the baseline that is enforced by 
mass conservation. For each pair of pulse distances and am- 
plitudes there is a well-defined velocity that is smaller than the 
speed of sound in such membranes. Further, the model implies 
that there is a minimum distance between pulses or a maximum 
frequency in a pulse train. Here, we have provided recordings 
from locusts (locusta migratoria) that display periodic pulses 
with long intervals, but with significant likelihood also min- 
imum distance pulse doublets. In contrast to other pulses in 
this nerve, the doublet spacing in time is unaffected by changes 
of the angle of extension of the locust leg, and no pulses of 
closer distance were found in the individual neurons. Since the 
femoral nerve is only about 2 cm long and the pulse propagation 
velocity is about 1 m/s (see results section) one can conclude 
that action potentials with a distance larger than 20 ms cannot 
be found simultaneously in a neuron. Therefore, they cannot 
influence each other and must be considered uncorrected. In 
the femoral nerve we found to classes of pulses. One class 
with large time intervals indicating separations too large for 
temporal coexistence, and pulse doublets with spatial distances 
smaller than the overall length of the nerves. Such pulses can 
influence each other, and we consider them here as an approxi- 
mation to a periodic solution for a nerve that has a length that is 
not larger than two pulse distances. The distance of the pulses 
in the femoral nerve is about 1 1 ms corresponding to about 1 cm 
in space. The individual pulse width is 1-2 mm (or 1-2 ms). 
Thus, the minimum pulse distance in the doublets is about 5-10 
pulse widths. The minimum distance in the electromechanical 
model is about 8 pulse widths. It should be stressed that the the- 
ory does not allow for stable doublet solutions. However, we 
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consider the doublets as an approximation to such solutions for 
the case of a nerve with a length only slightly larger than 2-3 
pulse intervals. It should be added here that in our electrome- 
chanical description a simple linear relation between density 
and voltage has been assumed. Therefore, the integral of volt- 
age changes over time (or space) should depend in a simple 
manner on the integrated density changes. 

Another feature of our theory is that it attributes a meaning 
to myelination. We have chosen parameters leading to pulse 
velocities of approximately 100 m/s, comparable to the pulse 
velocity of myelinated nerves and much larger than those found 
in the locust femoral nerve. We have noted before that our 
electromechanical theory assumes longitudinal compression 
Q), and we do not consider transverse excitations normal to 
the nerve membrane. However, it is well-known that longitu- 
dinal wave propagation is significantly faster than transverse 
propagation. E.g., the speed of sound in bulk water is 1500 m/s 
while a surface wave displays a typical velocity of about 1 
m/s. Therefore, it is to be expected that surface excitations in 
nerves involving changes in membrane curvature are consid- 
erably slower than longitudinal pulses. Thus, if one prevents 
out-of-plane waves by spatial confinement, one expects faster 
pulse propagation. This is the role of myelinization in the 
soliton theory. It has already been noted by Abbott et al. 
1 10 1 that the magnitude of the reversible heat changes suggest 
that the myelinated regions of the nerve play an active role 
in the pulse propagation process. This is not the case in the 
Hodgkin-Huxley model in which myelinization plays the role 
of an electrical insulator preventing an active amplification of 
the nerve signal. Therefore, they assigned a special role to the 
'nodes of Ranvier' which are the only parts of the myelinated 
nerve that could amplify the nerve signal. In the soliton model, 
the 'nodes of Ranvier' do not play an active role different from 
the rest of the nerve. The construction of a theory for transverse 
pulse propagation is a far more demanding task. However, such 
a model would offer a better description of non-myelinated 
nerves. Thus, although the present results are suggestive, their 
relevance to the non-myelinated locust femoral nerve has not 
been demonstrated. 

A further notable feature of the electromechanical soliton 
theory is that it does not require ion channel proteins. Such 
proteins are, however, indispensable in the Hodgkin-Huxley 
model. In the present model the most important requirement 
for obtaining localized pulses is the presence of a cooperative 
melting transition in the biomembrane. This transition is not 
only responsible for the reversible heat but also for the localiza- 
tion of the pulse. Interestingly, we and others have shown that 
in this transition the lipid membrane displays quantized ion 
conduction events that are indistinguishable from the channel 
events typically attributed to protein ion channels, even in 
the complete absence of proteins. In particular, one finds a 
straight-forward explanation of these quantized current events 
in inevitable thermodynamic fluctuations. (For a review, see 
[37]). Thus, during the pulse propagation, one expects that 



the lipid membrane should become permeable for ions and 
that quantized conduction events should occur. Similarly, any 
agency that inhibits pulse formation in the soliton model would 
also inhibit quantized currents since the two phenomena are 
unavoidably coupled thermodynamically. 

Summary 

We have shown here that the soliton theory for electromechani- 
cal action potentials in nerves possesses the following features: 

• It generates localized voltage pulses while displaying a 
reversible heat and mechanical changes including short- 
ening and swelling 

• When the overall length of the nerve is constant, one ob- 
tains periodic solutions. 

• One finds an undershoot or hyperpolarization. 

• The periodic pulses display a minimum distance of about 
5-10 pulse widths as a consequence of mass conserva- 
tion. 

• This pulse separation is approximately the same as that 
found in the femoral nerve of the locust. 

• During the pulse one finds quantized ion conduction 
events through the membrane resembling those usually 
attributed to ion channel proteins. 

All features of the nerve pulses and the generation process dis- 
play unavoidable thermodynamic couplings with a predictable 
influence of changes in temperature, lateral pressure, pH, cal- 
cium concentration or anesthetics. Thus, we have shown that 
the soliton theory for nerves is able to explain most known fea- 
tures of nerves, the effect of anesthetics [23] and the emergence 
of ion channel phenomena from the thermodynamics of the 
membrane without employing molecular features of the mem- 
brane components. 
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